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Abstract 

Actomyosin networks are major structural components of the ceh. They provide mechanical 
integrity and allow dynamic remodeling of eukaryotic cells, self-organizing into the diverse patterns 
essential for development. We provide a theoretical framework to investigate the intricate interplay 
between local force generation, network connectivity and collective action of molecular motors. 
This framework is capable of accommodating both regular and heterogeneous pattern formation, 
arrested coarsening and macroscopic contraction in a unified manner. We model the actomyosin 
system as a motorized cat's cradle consisting of a crosslinked network of nonlinear elastic filaments 
subjected to spatially anti-correlated motor kicks acting on motorized (fibril) crosslinks. The phase 
diagram suggests there can be arrested phase separation which provides a natural explanation 
for the aggregation and coalescence of actomyosin condensates. Simulation studies confirm the 
theoretical picture that a nonequilibrium many-body system driven by correlated motor kicks can 
behave as if it were at an effective equilibrium, but with modified interactions that account for 
the correlation of the motor driven motions of the actively bonded nodes. Regular aster patterns 
are observed both in Brownian dynamics simulations at effective equilibrium and in the complete 
stochastic simulations. The results show that large-scale contraction requires correlated kicking. 
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I. INTRODUCTION 



The mechanical integrity of eukaryotic ceUs depends on their cytoskeleton. The cytoskele- 
ton is made up of a dense network of protein filaments spanning the cytoplasm. Cytoskeletal 
networks self-organize into highly dynamic and heterogeneous patterns from the interplay 
between active force generation by molecular motors and passive dissipation of energy in the 
crowded cellular interior.- Understanding the dynamics of such pattern formation remains 
a challenge to statistical mechanical theory. 

Actomyosin networks are the main components of the cellular contractile machinery es- 
sential for processes as diverse as cytokinesis and wound healing. Walking on the structural 
scaffold provided by an actin network, myosin-II motors themselves self-assemble into bipolar 
minifilaments that generate sustained sliding of neighboring actin filaments relative to each 
other. By carrying out this correlated motion the minifilaments reorganize the filamentous 
actin networks and generate tension ultimately powered by ATP hydrolysis. The forma- 
tion and coalescence of actomyosin aggregates to exert contractile forces are manifested in 
pulsed contractions of an actomyosin network that drive epithelial sheet deformation during 
morphogenesis. Such aggregates also are responsible for a multistage coarsening process 
that occurs in a bottom-up model system for contractility which has been reconstituted in 
vitro.- 

In reconstituted filament-motor assemblies, relatively regular patterns such as asters, 
in which stiff filaments or filament bundles radiate from a common center arise that 
resemble the mitotic spindles formed in dividing cells. A polarity sorting mechanism^ has 
been proposed to explain the observed pattern. On the other hand, irregular heterogeneous 
cluster structures have been seen both in the actomyosin networks of C. elegans embryos^ 
and in the minimal in vitro network model. ^ In the in vitro system, clusters of various 
sizes continuously change via fusion and rupture events but the distribution of cluster sizes 
remains steady. Motor-driven filament sliding and a well-tuned connectivity seem to play a 
key role in this nonequilibrium steady state. 

In the present work, by looking into the intricate interplay between local force generation, 
network connectivity and collective motor action, we seek to provide a theoretical frame- 
work that is capable of accommodating both regular and heterogeneous pattern formation, 
as well as arrested coarsening and large-scale contraction in a unified manner. Our micro- 
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scopic model is a motorized cat's cradle^i""— which consists of a crosshnked network of non- 
hnear elastic filaments where motors generate equal but oppositely directed kicks at motor- 
bonded node pairs. This model encodes two essential aspects of actomyosin self-organization: 
first that individual actin filaments have highly asymmetric responses to loading, resisting 
large tensile forces but easily buckling under compressive loads. As a consequence, sheared 
biopolymer networks exhibit negative normal stresses^^ comparable in magnitude to the 
shear stress. Secondly the model captures the fact that motor-induced node movements 
come in anti-correlated pairs owing to the bipolar minifilaments. This anti-correlation mim- 
ics the relative sliding of neighboring parallel filaments driven by the attached motors. Both 
aspects are crucial to capturing the formation of a disordered condensed state of actomyosin 
aggregates. 

Our coarse-grained model is completely microscopic but bears some resemblance to the 
macroscopic approach adopted by Levine and MacKintosh that introduces force dipoles into 
an elastic continuun>i^ii^ where motor unbinding kinetics leads to enhanced low-frequency 
stress fiuctuations. In the microscopic model the network connectivity and motor distri- 
bution over the bonds of the network are quenched once initially assigned, so that the 
nonequilibrium dynamics and structures predicted by our model arise solely from the intrin- 
sic activity of motors firmly built into the network driving correlated motions stochastically. 
This assumption is in line with the fact that the in vitro structures are irreversibly assem- 
bled because many protein factors found in vivo that allow fast pattern renewals are left 
out of the reconstitution, such as disassembly of contractile structures and transience of 
actin crosslinking proteins. The model highlights the key role in determining the course of 
structural development played by the motor susceptibility, a parameter characterizing how 
sensitively the motors respond to imposed forces. The coupling between motor kinetics and 
the structure leads to a double-way feedback: Motor action induces structural changes of 
the network and thus modifies the local mechanical environment of the motors, which in 
turn changes the load-dependent motor response. 

In the same spirit as our earlier work on the statistical mechanics of systems with uncor- 
related kicks on each node^^^ we adopt a master equation description but treat small-step 
spatially anti-correlated kicks in order to mimic contractile-ratchet-like^'^ incremental de- 
formations of actomyosin networks. We again obtain an effective Fokker-Planck equation in 
the small kick limit. But the correlations lead to local effective temperature Tgfr and diffu- 
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sion coefficients Dgff, which now both become tensors and also depend on the instantaneous 
local network structure. More interestingly, anti-correlated kicking leads to a modification 
of the bare interaction. Essentially new forces come into play through the action of the 
motors. This motor-induced force depends linearly on the motor activity to quadratic order 
in kick step size and decays in space as the inverse distance in three dimensions, resembling 
a logarithmically growing potential. By treating the motor-bonded node pairs as "functional 
units" and deriving a pair-level steady-state solution of the effective Fokker-Planck equa- 
tion, we demonstrate that the motor-driven anti-correlated movements of actively-bonded 
nodes give an additional effective pair potential that exhibits a strong short-range attrac- 
tion regardless of motor susceptibility. The original interaction is also still present but is 
an enhanced or weakened long-range attraction. For motors with negative susceptibility at 
sufficiently high activity, the dominant interaction is actually a long-range repulsion. The 
relative contributions of the various terms yield a diverse range of steady-state structures. 
This decoupling scheme in the mean-field spirit also allows us to perform a self-consistent 
calculation to evaluate quasi-thermodynamic phase diagrams. A non-monotonic dependence 
of the pressure (or tension) upon the node density (under susceptible motor kicks) indicates 
the possibility of phase separation. 

To test the validity of the analytical approach, we compare the steady-state structural fea- 
tures found by Brownian dynamics simulations using the effective temperature and modified 
potential (both a tensor-parameter formalism and a pair-level scalar-parameter formalism) 
obtained from the steady-state solutions, with the structural features observed in a dynamic 
Monte Carlo simulation that is fully consistent with the master equation at thermal tem- 
perature and with bare interactions. We find good quantitative agreement suggesting that 
a non-equilibrium system driven by small-step correlated motor kicks can be thought of as 
being at an effective equilibrium with modified interactions. 

The existence of an effective short-range attraction combined with the predicted tendency 
for phase separation suggests that the formation of steady heterogeneous cluster structures 
is an example of arrested phase separation.— A force-percolating network consisting of 
nonlinear elastic fibers attains rigidity when local collapse induced by the motor-driven 
short-range attraction balances the concomitant neighboring bond stretching. This initially 
homogeneous network then develops into dense clumps connected by highly stretched bonds, 
and simultaneously, compact aggregates phase separate from node-poor regions. The coars- 
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ening process stops once a global balance is achieved, and the pertinent dynamic process 
involves phase separation followed by arrest due to bond constraints. Our model naturally 
explains the formation of aster patterns through the notion of an effective repulsion. 



II. THEORY 

In our earlier work^^ we showed how an effective temperature describes steady-state 
fluctuations and responses of a model cytoskeleton, treated as an amorphous network of 
crosslinked nonlinear-elastic filaments, driven by uncorrelated motor kicking events. Here 
we consider anti- correlated kicks. As sketched in Fig. [U each motor (myosin minifila- 
ment) generates a pair of equal but oppositely directed displacements (red arrows) at the 
motor-connected crosslinks/nodes (purple spheres). These anti-correlated kicks mimic the 
contractile-rachet-like incremental movements due to myosin-driven relative sliding of neigh- 
boring actin filaments. We point out that myosin motors do not explicitly enter our model; 
instead, they are exemplified only through the anti-correlated kicks. The cartoon in Fig. [1] 
illustrates how these kicks are generated: A myosin minifilament attaches to two otherwise 
unconnected actin filaments and pulls the node on either filament toward each other. We 
then assume, for simplicity, that an implicitly-motor- attached filament/bond connects the 
node pair, neglecting the detailed architecture of the motor-filament composite (as seen in 
the zoom-out view in Fig. [T]top image). Anti-correlated kicks then act on the nodes at the 
two ends of such a (motor-attached) active bond. 

The asymmetric load response of individual actin filaments is encoded via a nonlinear- 
elastic interaction between the bonded nodes, defined by the pair interaction potential 
/3[/(r) = B(r — Le)Pj{r — L^Y /2. Here B(-) is the Heaviside step function and /37 gives the 
effective stretching stiffness of the filaments with (5 = l/ksT. An energy cost arises only 
when the contour length r of a bond exceeds its relaxed length Lg. We call this interaction 
a "cat's cradle" interaction.-"— For simplicity, we assume that all the bonds, no matter 
motor- attached or not, have the same relaxed length. Since the motor-driven forces exceed 
by far the piconewton-threshold of affordable compressive loads, the induced buckling within 
a percolating actin network gives rise to a tensegrity structure composed of collapsed and 
stretched elements (illustrated in Fig.[2j): In a permanently crosslinked network of filaments, 
such as the in vitro reconstituted networks, active sliding of filaments is constrained by pas- 
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FIG. 1. Schematic illustration of the spatially anti-correlated kicks acting on motor-bonded node 
pairs. (Central image) A bipolar myosin minifilament pulls in slack locally, generating a pair 
of equal size (/) but oppositely directed displacements (red arrows) at the motor-bonded nodes 
(purple spheres) along their line of centers, where r^j is a unit vector pointing from node i to 
node j. Upon zoom-out, this represents a typical functional unit (marked by a dashed circle in 
the top image) that generates incremental contractions within a crosslinked filamentous network. 
An enlarged view of the actin filament (bottom image) reveals its segmented structure. The size 
/ of the subunits determines the magnitude of the relative node displacements due to contraction 
events of myosin sliding. / is thus taken to be the step size of anti-correlated kicks in our model. 



sive crosslinking, in other words, local filament or bundle contraction is balanced by the 
stretching of neighboring filaments. An initially homogeneous network typically then devel- 
ops into dense floppy clumps (concentrated short green wiggly lines) connected by highly 
stretched filaments (long red straight lines). This phenomenon leads to the formation of 
disordered actomyosin condensates^ and can lead to active contractility.— 
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FIG. 2. Cartoon of the tensegrity structure composed of collapsed and stretched elements. In 
a crosslinked network of filaments, active filament sliding is stabilized by passive crosslinking. A 
tensegrity structure is formed once a global balance between local contraction and neighboring bond 
stretching is achieved. An initially homogeneous network then develops into dense floppy clumps 
(concentrated green wiggly lines) connected by highly stretched filaments (long red straight lines) . 



A. Quadratic expansion of the master equation: effective equilibrium with modi- 
fied potential 

1. Fokker- Planck (FP) /Smoluchowski equation for Brownian particles 

Consider a collection of Brownian particles (i.e. nodes of our model network) labeled 
with positional coordinates fi {i = 1, • • • ,A^). For infinitesimal step Brownian motion, 
the configurational probability density ^^({r^};t) is locally conserved and thus satisfies a 
continuity equation /dt = — Vf. • J^, where Ji is the probability current density along 
the coordinate of the ith particle. 

For Brownian particles the probability current density is linearly related to the deviation 
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of the configurational probability density from its equilibrium value 

J^ = -Y. Dliin}) ■ (V,-* + ^mv,u) . (1) 

3 

Here U{{ri}) is the thermal equilibrium potential of mean force of the system and (5 = 
{kBT)~^ . The equilibrium distribution ^^eq({^i}) is related to U ({r^}) according to ^^eq({^i}) oc 
exp[— /3[/({r^})]. The diffusion coefficients D^-^ are functions of the system's configu- 
ration and satisfy a generalized Einstein relation with the drag coefficients reading 
D^ij = kBT{(^)^j^^ where D^j and C-j are 3x3 matrices for each (i^j) pair. 

When divided into self-diffusion and coupled-diffusion parts, the FP equation /dt = 
LP ^ becomes 

FP 

i 

+ E E • ■ (^^-^ + /5*v,f/). (2) 

i j^i 

Note that the double gradient operation acts explicitly as 

V, • Dl ■ V,* = (V, • Dl) ■ V,* + Dl : V,V,*, 
Vi ■ Dl ■ W^VjU) = {Vi ■ Dl) ■ W^VjU) + Dl : {/3Vi^VjU + /3^ViVjU). (3) 

2. Master equation for motor- driven processes: anti- correlated kicks 

To mimic the motor-driven filament sliding in actomyosin networks, we describe the 
motors as generating anti-correlated kicks on pairs of crosslinks that pull in slack locally 
(Fig. [1] middle). Since the linear size of the myosin minifilaments is small compared to the 
mean separation between the crosslinks, the anti-correlated moves can be treated as being 
along the lines of centers. In view of the segmented structure of the act in filaments, which 
consist of periodically arranged subunits of linear size / (Fig. [1] bottom), we therefore assume 
a fixed kick step size /. / indicates the amplitude of relative node displacements due to a 
typical contraction event. Thus an anti-correlated kick pair acting on nodes i and j can be 
represented by a pair of displacements along the line of centers (/^j, Iji) = /(r^j, —fij)^ where 
fij is a unit vector pointing from node i to node j. These anti-correlated kick pairs with 
equal size automatically satisfy momentum conservation on the macroscopic scale. Yet if we 
include explicitly the aqueous environment in which the cytoskeletal network is immersed, 
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hydro dynamic interactions^^ between the nodes via the solvent should be taken into account. 
These interactions might modify the current simplified picture, and counteract any motor- 
induced force imbalance on individual nodes, thus validating momentum conservation on 
the microscopic scale as well. 

The dynamical evolution of the many-particle configuration {fj} due to these motor- 
driven events can be described by a master equation /dt = I/^e^ with 

L,,^{{fi}-t) = J Yidr\ [K{{r'^} ^ {{?,}■ t) - K{{n} ^ {r\})<^{{n}-t)\ , (4) 

where the integral kernel K{{r\} {^i}) encodes the probability of transitions between 
different node configurations. Our earlier description^^^^^ of the motor kicking rate, fc, still 
applies to current case for correlated kicks, i.e., 

k = /^[e(A[/) exp(-5^/3A[/) + e(-A[/) exp(-5rf/3A[/)], (5) 

where k is the basal kicking rate and Su{sd) denotes motor susceptibility to energetically 
uphill (downhill) moves, except that now the free energy change AC/ is due to pairs of 
displacements. Explicitly we write 

x^6{fi-Pi-kj)6{fj-Pj + kj)w U{'" ,r'i,'" ,r'j,'")-U{'" ,ri,'" ,rj,'") 
-5{ri - r'i + lij)5{rj - r' j - lij)w [/(••• ,fir" .rj," -) - U{" - ,r'i, • • • ,r'j, • • 

The factor 1/2 avoids double counting in the summation over all pairs. The quantity Qj, 
much like an element of a contact map in description of protein structures, defines whether 
the node pair {i^j) is connected by an active bond and thus subject to anti-correlated 
displacements (/^j, —lij)'- Cij = Cji = 1 for motor-bonded pairs while Cij = Cji = for non- 
bonded pairs. Our description of the rates gives w[Ui — Uf] = Q{Uf — Ui) exp[—Su/3{Uf — 
Ui)] + e{Ui - Uf) exp[-Sd/3{Uf - Ui)]. 

Assuming symmetric motor susceptibility, i.e. Su = = 5, one finds more simply 

i 3 

^ |e-«/3K.f,)-t/(ri-r,,/,+r.,)] . . ^ ^ - 4, ... , r', = r, + k,, . . . }; t) 
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We assume that kicks on different pairs of nodes at any given time are uncorrelated. The 
rates of possible kicking events depend on the instantaneous node configuration reflecting 
an assumed Markovian character of the motor dynamics. There is no angular average due to 
the deflniteness of kicking directions for a given conflguration. Note that the motor power 
strokes and thus the kick steps are discrete occurring in a stochastic fashion. The correlated 
motions pull in slack locally while simultaneously pulling taut neighboring fllaments until a 
global balance is reached or a macroscopic collapse occurs, depending on whether the motors 
are downhill prone (with a large positive s) or load-resisting (with a small or negative 5), 
respectively. 

Quadratic expansion of Eq. [6] in kick step size / leads to 



"1. . 1. . 

+s {fijhj ■■ ViVi/SU + fijVij : VjVj/SU - 2fijfij : ViVj/3U) ^ 



Notice that '^i'^j Cijrijfij : VjVj = '^i^j Cijfijfij : VjVj, the above expression can be 
rewritten as 



EE 



•(7) 



The deflnitions in Eq. [3] allow us to express I/^e^ form of the divergence of a flux 

plus some extra terms which modify the bare interactions as shown below. 



3. Generalized FP equation for motorized systems: effective temperature and modified potential 

By combining the pure Brownian dynamics (Eq. [2j) with the nonequilibrium dynamics 
due to correlated motor-driven processes up to 0(/^) (Eq. [Tj), we obtain an effective FP 
equation with effective tensor parameters [D^^ and (3^^ and sitewise modifled potentials 
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mod ) 



eff 



The effective diffusion constants read 



(8) 



(9) 



The effective temperatures are given by 



f3% = HDl 



SKjI Cij Tij Tij 



1 



(10) 



The modified potentials U^J^ = U + A*'^-'^ involve the (additive) modifications A that satisfy 



-V,;A^ = 



Cijfijfij^P 



.1 



VjA^ = Vi- {-KfCijTijrij) • {Dl - sKfCijfijrij)^ 



(11) 



We assume that the active connectivity described by {Cij} is quenched once initiaUy 
assigned and use the identity that Vj • {TijTij) = — • (rijTij) = rij{d — l)/rij^ where d is 



the spatial dimension and r 



r,- — r 



to obtain exphcit expressions. 



Several nontrivial features can be read off from the above expressions (Eqs. IS HTTI) : 
(1) There is a key difference from the simple situation for uncorrelated isotropic kicks.— 
In that case an average over kicking directions yields uniform scalars Teff and L^eff which 
are fully determined by the motor properties (motor susceptibility s and activity defined as 
A = tvl'^/Do) regardless of the specific structure of the system. Under correlated kicks along 
the lines of centers, however, both Teff and Deff become local tensors. These tensors depend 
on local network structure (relative position of bonded neighbors r^j) and motor distribution 
(quenched active connectivity defined by Cij) about the central node i. The tensorial na- 
ture of the effective diffusion coefficients (Deff) and mobility (/x^f^ = D^^ • /3^f^) leads to the 
diffusive flux not aligning with the density gradient or the drift flux. Motor-induced modi- 
fications to the forces (— Vit/mod) and to the transport coefficients (I?efr and jn^^) only have 
longitudinal components since fijfij essentially serves as a longitudinal projection operator. 
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(2) Anti-correlation of the myosin-generated kicks at the two ends of the actin filaments 
provides a microscopic basis for the anisotropy of actin diffusivity: The anti-correlated 
kicks enhance self-diffusion of individual molecules, whereas they slow the relative coupled 
diffusion of motor-bonded molecules (Eq. [9j). Recent experiments^S indeed have reported 
that motor activity can give rise to cytoplasmic motion that has the appearance of diffusion 
but is significantly enhanced in its magnitude. Such "active" cytoplasmic diffusion could 
enable rapid intracellular transport of matter and information. 

(3) To quadratic order in kick step size /, it is easy to show: 

(a) Both Vi (AY/c^r) and Vj (A-^/fc^T) are proportional to k.I'^/Dq. This indicates that 
the motor-induced forces depend linearly on the strength of the kicking noise relative to the 
thermal noise in small kick limit. 

(b) The inverse effective temperature tensor becomes 



= /3 



+ 0{l 



4\ 



Anti-correlated kick pairs give rise to effective temperature tensors that depend on local 
active connectivity, implying that the speed of heat flow in motor-kicking directions [Cij ^ 
thus Teff 7^ T) is different from energy flux along motor- free directions {Cij = thus Teff = T). 
Along individual active bonds, as for the uncorrelated kick case, we have Tgff <Tif5>l/2 
whereas Teff > T if 5 < 1/2, as well as the same detailed balance condition that if 5 = 1/2 
then Teff = T. 

Therefore the nonequilibrium system driven by small-step anti-correlated motor kicks 
may be described by an effective FP (Smoluchowski) equation at local effective temperature 
and with modified interaction potential. 



B. Pair- level steady-state solution 

We may ask whether it is still possible to obtain an explicit (but perhaps approximate) 
steady-state solution to an effective Fokker-Planck equation with tensor transport coefficients 
and modified interaction potential. Existence of such a solution makes possible the mapping 
of the system to an equilibrium system even for this nonequilibrium situation with spatially 
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correlated motorized events. This mapping thus allows the study of rheological properties 
within a quasi-equilibrium framework. The main difficulty in making this mapping arises 
from the tensorial nature of the effective temperature which depends on local structures. 

To approximate the A^-body solution, we employ a decoupling scheme that reduces the 
problem to finding the steady states for the diffusion of "functional units". We will first 
study the simplest case for a single pair of motor-bonded nodes without hydrodynamic 
interactions, i.e., neglecting the D^j-related terms. A straightforward calculation leads to 
an explicit steady-state solution to this two-body problem. 

Consider a pair of nodes located at ri and r2- For purely Brownian motion, the 2-body 
configurational probability density ^^(ri,r2;^) evolves according to the bare Fokker-Planck 
equation /dt = ^^p^^ with 

Here the mobility tensor is related to the diffusion tensor simply by a multiplication of the 
inverse scalar temperature as jj.^- = /3D^-. The steady-state solution thus obeys the usual 
Boltzmann law ^^^^ oc exp[— /3[/(ri2)] where the interaction potential U only depends on the 
node separation ru = \ fi — in the absence of external fields. 

The dynamic evolution of a pair due to anti-correlated motor kicks is described by a 
master equation /dt = I/^e^ with 

4,^^(fi,f2;t) = ^Cule-'^^""^''^''^^^^^^ = n - h2.r2 = r2 + h2;t) 

By combining the pure Brownian dynamics (Eq. [T2l ) and the motor-driven dynamics 
(Eq. [13]) we obtain in the small kick limit an effective Fokker-Planck equation characterized 
by effective diffusion constants, tensor effective temperature as well as a modified potential 

i=l,2 

In this two-body case, the effective diffusion tensors, up to are given by 

Df, = Dol + Ci2^«:/'fi2ri2 = DqI + C2i\nl''r2ir2i = D'g. (15) 
The effective mobility tensors, up to OQ?), are given by 

ixf^ = 13 (Dol + Ci2snfh2h2) = /? {Dol + Csis/^/^fsirsi) = ^2! • (16) 
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The tensor effective temperatures thus read 
of = f3f = P 



The modifications to the bare forces are central forces that decay in space as 1/r: 

- [-ViA] , 



ViA = Ci2^— ri2 , -V2A = Cai^— rai 



suggesting the following form of the modified potential 

U^odiru) = U{ri2) + Ci2—kBT\a{ri2) + const. 



(17) 



For s = 1/2, the effective temperature becomes uniform scalar again f3^^ = /31. The steady- 
state pair solution thus has precisely the form ^^ssi^u) oc exp[— /9[/inod(?'i2)] where t/mod is 
given by Eq. [T3 

For s 7^ 1/2, we have 

1 . 



/3'^-[-Vif/„od] = /9 



-. ^ , 1\ ■ 

^ + ' " " 2 ] Wo'''''' 
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f3Do ri2 



ri2 



1 + C12 u 



1\ 



2 Dn 



kP 1 

[-Vif/]+Ci27^— ri2 + 0(/4) 

^0 ^ 12 



and 



(18) 



(19) 



/3^" ■ [-V2f/mod] = -/3^" • [-Vif/^od] . 

Here we have used the identity fi2(fi2 • Vit/) = Vif/, since ViU is parallel to fi2 in the 
two-body situation. We see therefore that we still arrive at a steady-state solution given by 



"^ss oc exp < - j3 



1 + C12 s 



1\ 



2) Do 



U - Cu^^ Hrn) } = exp [-^esUes] , (20) 



where the scalar inverse effective temperature {/3es) and the effective interaction potential 
(t/eff) are given by 



and 



eff 



U + 



C12 (kP/Do) 



(21) 



(22) 



1 + ^2(5-1/2) (tiP/Do) 
Since on the pair level the total force is automatically along the line of centers, Teff 
becomes effectively a scalar. Thus at the pair level the steady states take a form consistent 
with our earlier result for uncorrelated isotropic kicks^i (note that here is no factor that 
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arises from the angular average). The modification to the bare interaction now, however, 
yields an additional central force which decays in space as 1/r. In the limit of high motor 
activity, i.e., A = nf/D^ > 1, Eq. [22] reduces to U^^ - U c:^ kBTln{ri2)/{s - 1/2). We 
see the motor-induced forces can be attractive or repulsive depending on whether the motor 
susceptibility s is larger or smaller than 1/2, respectively. 

The scaled effective pair potential given by ^esUeff = [1 + Ci2{s — l/2)A]+(7i2A ln(ri2) 
consists of two terms. The first term indicates that the motor action may enhance or weaken 
the long-range attraction, arising from the ordinary bond stretching in the model cytoskele- 
ton. The sign of this effect depends on the sign of (5 — 1/2). The second effect of the 
anti-correlated motors is an effective confinement potential that promotes further contrac- 
tion even in the buckling regime. This term is independent of motor susceptibility. Therefore 
by varying s we may distinguish the influences of either term. 

Fig. [3] displays the profile of the effective pair interaction [/eff (Eq. [22]) scaled by the 
effective temperature (Eq. [21]). As shown in panels (a) and (c), when the motor suscep- 
tibility s is varied the long-range interaction becomes modified but the effective attraction 
at small separations is not affected. In contrast, increasing the motor activity A not only 
influences interactions at large distances but also enhances the short-range attraction, as 
seen in panels (b) and (d). Notably, at sufficiently high activity, load-resisting motor kicks 
{s < 0) may yield a long-range repulsion, i.e. the slope of the ^^efft/eff curve becomes negative 
at large r (see two lowest curves in panels c and d). As a consequence, an energy barrier 
(indicated by the red arrow) appears at intermediate distances, indicating the tendency for 
node separation and thus bond stretching that ultimately leads to aster formation observed 
in simulations. 

C. Self-consistent phonon (SCP) calculation: possibility of phase separation 

The pair-level steady-state solution derived above allows us to obtain the effective pair 
potential and self consistently determine the Debye- Waller factor of the localized nodes 
using the self-consistent phonon (SCP) method. It is straightforward to derive the effective 
potential associated with a "cat's cradle" with excluded volume, i.e., a system where the 
nonlinear bonds (with an effective stiffness /3^) stretch elastically when their contour length 
exceeds the relaxed length Lg but that buckle upon shortening too much. We assume a 
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FIG. 3. Profile of the modified interaction given by the pair-level steady-state solution. We plot 
the effective interaction 17^^ (Eq. [2T|) scaled by effective temperature (Eq. [22]) for various motor 
activity (A) and susceptibility (s). At sufficiently high activity, load-resisting {s = —0.3) motors 
may yield a long-range effective repulsion, an energy barrier (indicated by a red arrow in panel d) 
thus appears at intermediate distances, suggesting the tendency for node separation and thus the 
bond stretching that occurs in aster formation, (a) and (c): s = 1,0.5,0, —0.2 and —0.3. (b) and 
(d): A = 0, 0.5, 1, 2 and 4. Common parameters are Lg = 1.2, /37 = 5, Pc = 0.4, and Pa = 1. 



negligibly small hard core size (a) for comparison with simulations. We assign two mean- 
field parameters to characterize the network architecture: (1) network connectivity, Pc, 
which denotes the fraction of nearest-neighbor pairs that are bonded by filaments; (2) motor 
concentration, P^, which indicates the fraction of active bonds that are attached by motors 
and thus induce anti-correlated kicks on the connected node pairs. 

By using the independent oscillator approximation which yields sitewise decoupling of 
the particles, the free energy (due to configurational degrees of freedom) is expressed as a 
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sum of the effective potentials between the interacting density clouds 

/JKffdn - ffla,) = -Iny dr,p,(r;)e-^^^(^-'-^~^), (23) 

which essentially averages the Mayer f-bond, exp [— (l/2)/3[/(r^ ~ ^j)]^ over the location of 
the fj particle with an assumed Gaussian density distribution about the fiducial position Tj 

pirl,ri)=(^Y\-'^^^^^^-^i^\ (24) 



where d is the spatial dimension of the system. 

Self consistency requires that the eflFective potential mimics the harmonic comparison 
potential. Thus the curvature of the effective potential must coincide with the phonon 
frequency or spring constant of the Einstein harmonic oscillators. This gives a coupled set 
of self-consistent equations for {a^} 

«^ = ^ E Tr[VVV,si\r! - rj\;a,)]. (25) 

j 

In the present work the force constants {ai} will be taken to be equal. This is an effective 
medium approximation. 

For a cat's cradle with excluded volume and subject to correlated motor kicks described 
by the pair-level interaction (Eqs. [2TH221) . the effective potential /3efrKfr can be expressed as 



e 



a 1 r 
^RJo 



X f{w), (26) 



where the function f{w) encodes the pair interaction depending on whether the nodes are 
bonded and whether the bond is motor-attached. f{w) can be written in terms of the 
Heaviside step function Q(x) as follows 

f{w) = (1 - P,)e{w - a) X 1 

+Pc{l - Pa) \0{w - a)e{Le -w)xl + e{w- L^) X e-i^-^^^-^-)' 

+P,Pa \Q{w - a)Q{L, -W)X ^ _ ^ ^-A/2g-i/?7[l+(.-l/2)A](»„-L.)2 

(27) 

Here the first line denotes the interaction between non-bonded pairs where only hard-core 
repulsion enters (a stands for the hard-core diameter); the second line is due to the interac- 
tion between passively-bonded pairs where elastic bond stretching takes place beyond L^] the 

17 



third line accounts for the interaction between actively- bonded pairs where anti-correlated 
motor kicks induce an effective attraction even in the buckling regime {a < w < L^.) and 
yield an effective bond stiffness that depends on motor activity and susceptibility. 




FIG. 4. Profile of the effective pair interaction (Eq. [26l) obtained by the self-consistent phonon 
calculation, (a) Susceptible motors (s = 1) with increasing activity A enhance the long-range 
attraction and strengthen the short-range effective attraction, (b) A zoom-in view of the small- 
R region close to the elasticity onset (dashed line) in panel (a) showing the absence of kink or 
inflection point in the potential profile, (c) Adamant motors {s = —0.3) weaken the long-range 
attraction. No stable a solution is found if the motor activity gets too high (A > 1). (d) Varying 
motor susceptibility does not affect the short-range effective attraction (inset), but increasingly 
susceptible motors (bottom to top) lead to a stronger long-range attraction. Common parameters 
are = 1.2, f3-f = 5, = 0.4, Pa = 1. 



Fig. m shows the profile of the effective potential /3efrV^fr(i?, cx) given by Eq. [26] at various 
values of motor activity (A) and susceptibility (s). Panel (a) illustrates the logarithmically 
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growing potential at small separations due to motor-induced effective attraction, and the 
quadratically increasing potential at large distances due to ordinary bond stretching. Under 
susceptible motor kicks (5 = 1), both the short-range and long-range attractions become 
enhanced as the motor activity increases. The dashed line marks the location where the 
elasticity of individual filaments sets in (i.e. R = Lg). A zoom-in view of the small- i? 
region in panel (a) close to the elasticity onset (shaded region) is presented in panel (b). 
This figure clearly shows that the effective average potential increases monotonically with 
increasing separation; there is no inflection point or kink in the potential profile which 
would kinetically slow binding. Thus no difficulty is expected for initial aggregation, as ev- 
idenced by simulations. For load-resisting motors (panel c), however, as the motor activity 
increases, long-range attraction apparently becomes weakened due to a higher Teff, while 
the susceptibility-independent short-range attraction still gets stronger. No finite a solution 
can be stabilized when the motor activity gets too high (A > 1), signaling the develop- 
ment of spatial heterogeneity. Inhomogeneous/Site-dependent a solutions should recover 
the energy barrier at high motor activity, as seen for the pair-level solution (Fig. 3c and d). 
We demonstrate in panel (d) how the effective potential changes with motor susceptibil- 
ity. As expected, varying s has little effect on the effective attraction at short distances 
(more clearly seen in the inset), yet increasingly susceptible motors greatly strengthen the 
long-range attraction. 

In our model, given the asymmetric load response of the nonlinear elastic bonds, the 
effective attraction between motor-bonded nodes may well drive phase separation of a force- 
percolating network into dense clumps against voids, i.e. local condensates of contracted 
bonds (analogous to the droplets of the condensed phase) connected by stretched bonds (as- 
sociated with surface tension of the droplets), as sketched in Fig. [2l Spontaneous formation 
of disordered aggregates has indeed been observed in reconstituted actomyosin systems.-*^ 
Computer simulations of the model network when subjected to correlated motor kicks also 
show that substantial spatial heterogeneity develops when we started from a uniform distri- 
bution of motors over an isotropic network (see Fig. [5^ for an example). 

One thermodynamic indicator of phase separation is a non-monotonic dependence of the 
pressure upon the density of the constituents. We thus examine how the pressure varies 
with the node density using the self-consistently determined phonon frequency. To avoid 
structural complexity we perform the calculations on regular lattices, yet we expect the same 
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qualitative behavior for random lattices where an isotropic radial distribution function for 
the fiducial configurations can be obtained from simulations. The pressure p for a simple 
cubic lattice with a lattice constant R and a total number of nodes at the ambient 
temperature T can be evaluated by numerically diflFerentiating the free energy F as given by 



where F = z/3effV^ff (i?, Q^; A, s) with coordination number z = 6 for a simple cubic lattice. 
The dependence of pressure upon motor properties arises from the effective interaction 
encoded through f{w) (Eq. [271) and from the self-consistent a solution. 

[Note that since we have ignored the infiuence of motor-induced effective interaction on 
the fiducial structures, the SCP calculation may overestimate the instability threshold of the 
homogeneous state in terms of motor susceptibility s. Nevertheless, given the perturbative 
nature of the pair solution, any quantitative deviation should be modest.] 

III. SIMULATIONS 

To test the validity of the idea of using local effective temperature along with a modi- 
fied potential for the motorized system, we have performed three types of simulations and 
compared the resulting steady-state structures. These simulation models are: 

(1) Brownian dynamics simulations using a ien^or-parameter Langevin equation consis- 
tent with the effective Fokker-Planck equation (Eq. [8j) 

(2) Brownian dynamics simulations at Tgff (Eq. [2T1 ) and with a modified potential [/gff 
(Eq. [22]) using a 5ca/ar-parameter Langevin equation based on the pair-level steady-state 
solution 

(3) Stochastic simulations incorporating anti-correlated kicks along individual active 
bonds as chemical reaction channels, at bath temperature T and with the bare interaction 
potential U{{fi}). This is of course the most realistic model. 

A. Simulation setup 

We investigate a model cytoskeleton that consists of nonlinear elastic filaments subjected 
to anti-correlated kicks on the motor-bonded node pairs. The bare interaction U{r) between 




(28) 



20 



bonded nearest-neighbor pairs is taken to be of the cat's cradle type defined earher, i.e. 
PU{r) = B(r — Le)/3j{r — Le)^/2. The assumed weakness of the excluded volume effect 
allows large-scale structural rearrangements to occur rather readily. 

We build the model network on a simple cubic lattice to avoid structural complexity, and 
connect the nearest-neighbor nodes with nonlinear elastic bonds at a given probability Pc. 
Note that despite the regular lattice structure, disorder is still inherent in the randomness of 
bond connectivity for a partially connected network {Pc < 1). In most of the cases that we 
will study, we assume = 1, i.e., all the bonds are motor- attached. The bond connectivity 
and motor distribution are quenched once initially assigned; there are no bond or motor 
rupture events. We choose the filament relaxed length Lg to be larger than the lattice 
spacing (set as the length unit in simulations), so that the initial homogeneous network 
is completely fioppy with no tense bonds at all. The system size is = 6^ and periodic 
boundary conditions are applied. 

Since we are interested in the steady-state behavior at an effective equilibrium, we have 
chosen a kick step size that is sufficiently small such that higher order contributions to 
the /-expansion are not significant. On the other hand, the kick size is large enough such 
that the motor-induced effective interaction out-competes the thermal randomization. What 
affects the dynamics and structural development is the dimensionless motor activity or Peclet 
number A = tvl'^/Do which describes the relative strength of the motor kicking noise with 
respect to the thermal noise. Since in Brownian dynamics formalisms the motor kicking 
rate ^ and the kick step size / always appear in combination as an appropriate kicking 
rate has been chosen such that A > 1 yet no instability occurs. [In stochastic simulations, 
however, a higher basal kicking rate tv yields a faster approach to the steady state without 
influencing the steady-state features, since the basal kicking rate does not affect the relative 
probability of different kicking events.] 

Brownian dynamics (BD) simulations^^ have been implemented via the position Langevin 
equation Ar^(t) = /x^^ • (— Vit/mod)^^ + where R^^{At) represents the random 

motion due to thermal noise. The tensor formalism is equivalent to the effective Fokker- 
Planck equation (Eq. [8j) yet without hydrodynamic interactions. The total mobility of node 
i is given by jul^^ = /3{Dol + stvl'^ Cijfijfij)^ and the modified interaction force (— V^t/mod) 
comprises the total mechanical force {—ViU) acting on node i and the effective interaction 
(— V^A) given by Eq. dH The movement due to motor-induced effective attraction thus 
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follows /Xgfj- (— V^A) = tvl'^ ^ij^ij/^ij to quadratic order in /. The scalar formalism, based 
on the pair-level steady-state solution (Eq. [20]), sums up the contribution from individual 
bonded neighbors i.e., /x^^ • (-V,[/^od) ^ /^eff ' ^lo^ = /^(^o + ajSd^)F'^ (r,,) + 
/^/^ Ej ^ij^ij/^ij^ where F'^^{rij) is the bare interaction force between node i and its bonded 
neighbor j. The stochastic dynamics governed by the full master equation (Eq. [6]) has been 
realized by implementing dynamic Monte Carlo simulations^^ that obey the model kinetic 
rate (Eq.Ej). 

For making comparisons between the various simulation schemes, we ensure that all runs 
utilize the same lattice structure, bond connectivity and motor properties. Converging 
steady-state behavior would then validate the picture of an effective equilibrium at the 
effective temperature and with the modified interaction as predicted by the expansion. 

B. Illustrations 

1. Validity of an effective equilibrium and arrested phase separation 

An explicit way to test the validity of picturing the non-equilibrium system driven by 
small-step motors as being at an effective equilibrium is to compare the steady-state char- 
acteristics resulting from the three types of simulations that we described above. 

At a modest kick step size (/ = 0.03), all three simulation schemes lead to quite similar 
steady-state behavior despite disparate dynamics toward the steady state. The structural 
characteristics include (1) the mean squared node displacement (MSD) with respect to 
the initial regular configuration and (2) the amplitude of the innermost peak of the pair 
distribution function (PDF) which reflects the strength of aggregation. The PDF is defined 
as PDF(r) = C5n{r^ r + Sr)/r'^6r where 5n(r, r + Sr) counts the number of particles within 
an interval 6r at a distance r from the central particle, and the numeric factor C takes care 
of normalization. Both of these measures are almost identical for all the three schemes; they 
saturate to the same steady-state plateau value (see Fig. [St, d). The PDF profile in the 
main panel of Fig. [5]1 is obtained by averaging over a wide steady-state time window. Note 
that the compact aggregation, reflected in the steep rise of the PDF in the vicinity of the 
central node (i.e. the dominant peak at the minimum separation), results from the motor- 
induced short-range attraction in addition to the absence of excluded volume. The inset 
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shows the time evolution of the aggregation strength. The shght disparity in the potential 
energy (Fig. [5^) and the fraction of taut bonds (Fig. [Sh) seen in each simulation might arise 
from the perturbative nature of the expansion as well as from the difference in dynamics. As 
the kick step size increases we would expect larger deviations. The close resemblance of the 
steady-state node configurations and bond structures between different schemes (Fig. [5^) 
lends explicit support to the equivalence of the three schemes in the small kick limit, thus 
validating the picture of an effective equilibrium with modified interactions. 

To demonstrate the arrested phase separation, as anticipated from our theory, the network 
connectivity has been chosen to be sufficient for tension percolation yet moderate enough 
to allow considerable local force asymmetry (Pc = 0.4, i.e. z = 2-3). The initial network 
is completely fioppy (in green) with no tense bonds (in red) at all and the nodes sit on 
a simple cubic lattice (Fig. [5^ left panel). As shown in the right panels of Fig. [5^, under 
anti-correlated susceptible motor kicks, the active nodes (those with motor-attached bonds; 
shown as red spheres) begin to aggregate and tend to separate from the passive nodes (those 
with no motor-attached bonds; shown as blue spheres). The corresponding network structure 
exhibits clumps of fioppy bonds (concentrated short green lines) connected by tense bonds 
(long red lines). The overall rigidity of the structure is protected by susceptible motors 
which tune the balance between local bond contraction and neighboring bond stretching 
such that energetically unfavorable tense states are avoided. 

The arrested phase separation builds up as follows: Susceptible motors enhance the po- 
tential gradient via an (^-dependent) effective temperature, since Teff < T along the kicking 
direction if 5 > 1/2. This leads to an enhanced long-range attraction and thus strengthens 
the initial trend of aggregation among the motor-bonded nodes. The motor-induced short- 
range attraction then efficiently makes the aggregates become compact, yielding a phase 
separation of an initially homogeneous structure into node-rich and node-poor regions. In 
analogy to the nucleation of liquid drops within an initially homogeneous gas, where large 
surface tension serves as the driving force to form a bulk condensed phase, in our picture, 
local aggregates correspond to the droplets of the condensed phase whereas the stretched 
bonds connecting the aggregates contribute to the surface tension. Therefore, an ensuing 
coarsening process serves to reduce the surface area via coalescence of local aggregates into 
larger condensates. Once a balance between local bond contraction and neighboring bond 
stretching (given a force-percolating network structure) is achieved, the coarsening pro- 
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FIG. 5. Testing the validity of the effective equihbrium approximation: a comparison of three 
simulation schemes. Statistical characteristics and steady-state structures for a partially and ran- 
domly connected (Pc = 0.4) network built on a simple cubic lattice driven by small-step (/ = 0.03) 
susceptible (s = 1) motors are shown, (a) The potential energy; (b) the fraction of taut bonds; 
(c) the mean square node displacement; (d) main: the pair distribution function (PDF) averaged 
over a wide steady-state time window; inset: the aggregation strength, which is the height of the 
innermost peak of the PDF, versus simulation time; (e) initial (left) and steady-state (right) node 
configurations (upper row) and corresponding bond structures (lower row). The parameters chosen 
for illustration are = 1.2, /37 = 5, — 1? and = 1. 



cess stops and the structure does not evolve any further (as reflected in the plateau of the 
aggregation strength shown in Fig.[5ll inset). The system ends up with an arrested heteroge- 
neous structure with compact aggregates/dense clumps coexisting with voids/dilute regions 
(Fig. [5^ right panels); only moderate fluctuations about the arrested structures have been 
observed. Therefore, local force asymmetry is necessary for the initiation of phase separation 
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while force percolation is essential for achieving global balance and thus the eventual arrest. 




FIG. 6. The dependence of network tenseness and structure on motor concentration (P^) obtained 
by Monte Carlo simulations. The parameters were chosen such that the system is in the regime of 
arrested phase separation, (a) The fraction of taut bonds decreases as Pa increases. A kink located 
around Pa = 0.7 separates two descending branches: (I) Pa = 0.5-0.7 and (II) Pa = 0.8-1. (b) 
The aggregation strength exhibits a sharp peak at Pa — 0.7. The error bars in (a) and (b) depict 
standard deviations from averages over a steady-state time window of 4 x 10^ Monte Carlo steps, 
(c) Bond structures and corresponding node configurations at various Pa values are shown, from top 
to bottom Pa = 0.2,0.5,0.6,0.7 and 1. The arrow indicates the bond structure with the strongest 
aggregation. The remaining simulation parameters are Lg = 1.2, /37 = 5,Pc = 0.4,/ = 0.05,5 = 1, 
and n — 0.1. 

In Fig. [6] we demonstrate how the concentration Pa (and thus the cooperativity) of sus- 
ceptible motors (5 = 1) affects the tenseness and structural organization of the network. 
Parameters are chosen such that the system is in the regime of arrested phase separation at 
sufficiently high P^. Shown in panel (a), the overall trend of a decrease in the fraction of taut 
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bonds as Pa is raised is apparent. This results from increasing cooperativity in local force 
generation. Particularly noteworthy, however, is the presence of a kink near Pa = 0.7, which 
separates two descending branches (I: Pa = 0.5-0.7 and II: Pa = 0.8-1 as marked in the 
figure). In stark contrast, the aggregation strength (i.e. the amplitude of the first peak of 
the radial distribution function) exhibits a highly non-monotonic dependence on P^, which 
is sharply peaked at Pa = 0.6 (see panel b). The total energy (not shown) shows a similar 
trend as Pa varies. To understand these trends, we visualize the steady-state bond struc- 
tures and corresponding node configurations in panel (c). These simulated configurations 
reveal two distinct regimes characterized by different "strategies" to reduce the number of 
taut bonds: At intermediate motor concentrations (corresponding to branch I), as Pa rises, 
an increasing number of fioppy bonds are formed at the aggregation centers, at the cost of 
fewer but even more strongly stretched inter-clump filaments due to more compact node 
aggregation, as most clearly seen for Pa = 0.6 in panel (c), where the corresponding aggre- 
gation strength of nodes reaches its peak value (indicated by arrow). At sufficiently high 
motor concentrations (branch II), however, cooperative motor processes tend to minimize 
the overall tenseness by buckling as many bonds as possible at the same time avoiding too 
strong stretching. The resulting structure is thus rich in moderately buckled bonds and con- 
sists of large fioppy clumps connected by sparse tense bonds (see Pa = I case in panel c). 
The aggregation strength decreases with increasing Pa in this highly cooperative regime. 
Such an interesting dependence of structural development on motor concentration vividly 
demonstrates the intricate interplay of local force generation and collective motor action. 

2. Effective attraction 

To demonstrate the effect of motor-induced short-range attraction, we choose a network 
connectivity that is below the percolation threshold {Pc = 0.2). In the absence of global 
force percolation, when driven by spatially anti-correlated motors kicks, the initially homo- 
geneous and entirely fioppy network (Fig. [Tt upper left panel) develops into isolated fioppy 
clumps (Fig. [Tb lower panels), regardless of motor susceptibility. Under susceptible motor 
kicks (5 = 1), it becomes evident that further contractions of the already buckled bonds 
(with typical spots marked by circles in Fig. [7b lower left panels) arise solely from the ef- 
fective short-range attraction due to motor processes. On the other hand when driven by 
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FIG. 7. An illustration of the motor-induced effective attraction for a non-percolating network 
(Pc — 0.2) at various motor susceptibilities, (a) The potential energy; (b) the aggregation strength; 
(c) initial (upper left) and later node configurations and bond structures for a control run with 
pure thermal motion (upper right), and for motorized systems with 5 = 1 (lower left), 5 = 
(lower middle) and s = —0.5 (lower right). Despite having different dynamics, similar steady- 
state structures with isolated floppy clumps are reached in each case, regardless of the motor 
susceptibility. The remaining simulation parameters are Lg = 1.2, /37 = 5,Pa = = 0.03, and 

A^ = 0.1. 



load-resisting motors {s = —0.5), the anti-correlation in movements causes collapse of the 
transiently stretched bonds (Fig. [Tt lower right panels). Adamant motor kicks {s = 0) 
lead to a similar transient structure (with fewer tense bonds though) and the eventual 
collapse (Fig. [Tt lower middle panels). Therefore in a non-percolating network driven by 
anti-correlated kicks, despite the very different dynamics due to differing motor susceptibil- 
ity, similar steady state structures are reached. The complete collapse of all the individual 
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clumps is characterized by a vanishing total potential energy after the initial transients 
(Fig. [7K) and a significant aggregation strength that saturates to a steady-state plateau 
(Fig. [7h) once isolated condensates form. 

In contrast, the control run with pure thermal motion presents a considerable fraction of 
taut bonds (Fig. [Tb upper right panel) and thus maintains a finite potential energy (grey 
curve in Fig. [7K). The bond structure and node configuration remain largely homogeneous, 
exhibiting modest density fiuctuations and a vanishingly small aggregation strength (grey 
curve in Fig. [7h). 

Note that the illustrations for effective attraction shown here are obtained using fully 
stochastic simulations. The effective Brownian dynamics schemes give similar steady state 
structures composed of isolated fioppy clumps only for susceptible motor kicking. For load- 
resisting motors that give rise to an effective long-range repulsion, however, a distinct be- 
havior is seen (detailed below), highlighting the significance of correlation in motion for 
structural development. 

3. Effective repulsion 

Another interesting case arises when the motor susceptibility becomes negative. In this 
case a negative effective temperature yields an effective repulsion at sufficiently high motor 
activity (refer to Eqs. [2T14221) . As shown in Fig. [8](a), starting with an entirely fioppy network 
(upper row) , enhanced bond stretching coming from thermally induced fiuctuations in bond 
length giving local force asymmetries, eventually gives rise to highly tense and ordered aster 
patterns (lower row) at steady state. Aster formation occurs when the effective repulsion that 
promotes node separation and thus bond stretching dominates over the effective attraction 
that drives the opposite; removal of the short-range effective attraction does not affect aster 
formation, but indeed does disrupt efficient aggregation (see Fig. [TTb). Consistent with our 
earlier results for uncorrelated kicks,— the effective Brownian dynamics simulations also give 
sustained aster patterns. These patterns cannot collapse due to the absence of pairwise anti- 
correlation in motion; since at each move in Brownian dynamics, an individual node sees 
only an effective potential due to all its neighbors, the pairwise correlation is virtually lost. 
Such pairwise anti-correlation is crucial for active contractility as we showed elsewhere.— 
This anti-correlation, however, is captured by complete Monte Carlo simulations where anti- 
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correlated move pairs are treated as reaction channels and executed at each MC step. 



(a) (b) 

Pc=0.2 Pc=0.6 K=0.1 K=0.5 K=l 




FIG. 8. An illustration of the motor-induced effective repulsion caused by load-resisting {s = —0.5) 
motors for various network connectivities and motor kicking rates, (a) Initial (upper) and steady- 
state (lower) bond structures are shown at low (left: Pc = 0.2) and at high (right: Pc = 0.6) 
connectivity with = 1. (b) Node configurations (upper) and corresponding bond structures 
(lower) at various motor kicking rates (left to right: k = 0.1,0.5, and 1) with Pc = 0.4. The 
common set of simulation parameters are given by Lg = 1.2, /37 = 5, — I? and / = 0.05. 



Fig. [8] also highlights the architectural and dynamical ingredients required for the for- 
mation of connected tight asters that consist of tense bundles radiating from the common 
center. Panel (a) illustrates the necessity of a sufficient network connectivity for force trans- 
mission and bundle compaction; at low connectivity (Fig. left) only individually separate 
tense bundles are formed. Panel (b) depicts that a high motor kicking rate is needed to 
defeat thermal spreading and thus to facilitate filament or bundle alignment. 

Aster formation finds a natural explanation in our model when we use the notion of the 
effective long-range repulsion that we have derived. In Fig. [9] we show the mean-field indi- 
cations, obtained by SCP calculations (described in section II. C), of how the concentration 
(Pa) of load-resisting {s = —0.5) motors affects the long-range interactions. As clearly seen 
in panel (b), as Pa increases the long-range attraction due to bond stretching considerably 
weakens. Accordingly, the localization strength a of individual nodes (panel a) and the 
tension {—p) within the network (panel c) are both strongly suppressed. At P^ = 0.8, the 
profile of the effective potential becomes almost fiat at large distances (red arrow in panel 
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mean separation R mean separation R mean separation R 

FIG. 9. Mean-field predictions of the effect of the concentration (P^) of load-resisting {s = —0.5) 
motors on long-range interactions, (a) The localization strength a of individual nodes. The local- 
ization at large separation R is considerably suppressed as Pa increases, (b) The effective potential. 
Increasing Pa weakens the long-range attraction; the potential profile actually flattens out (red ar- 
row) at Pa = 0.8 indicating a vanishing restoring force, (c) The overall tension {—p) vanishes at 
large R (red arrow) for high Pa. This suggests the tendency for contraction is counterbalanced by 
a motor-induced long-range repulsion. The simulation parameters are Lg = 1.2, /37 = 5,Pc = 0.4 
and A = 1. 

b), indicating a vanishing restoring force. Consistently, a becomes vanishingly small at large 
R (panel a) and the overall tension decays to zero (red arrow in panel c), because the bond 
constraints are hardly felt when load-resisting motors counteract the tendency to contract. 
Yet higher motor concentration still enhances the effective short-range attraction (small- i? 
region in panel c) as expected. 

4. Effect of motor activity and susceptibility in phase separation: mean- field indications 

By performing the self-consistent phonon calculation, we find mean-field indications for 
the phase separation observed in the simulations and experiments. In particular, such calcu- 
lation allows us to examine how the changes in motor activity affect the tendency to phase 
separate which can be tested against simulations. We choose the length unit to be 100a 
(such that the excluded volume effect plays a negligible role in phase separation) and vary 
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the mean separation R between the nearest neighbors, i.e. the lattice spacing of the simple 
cubic lattice. 




FIG. 10. The effect of motor activity (A) and susceptibility {s) on phase separation. Shown are the 
calculated localization strength a (upper row) and the tension {—p) (bottom row) as a function of 
the mean separation R for (a) 5 = 1 with various motor activities, A = 0, 0.5, 1, 2, and 4 (bottom 
to top), and (b) A = 1 with various motor susceptibilities, s = —0.3, —0.2, 0, 0.5, and 1 (bottom to 
top). The remaining simulation parameters are Lg = 1.2, /37 = 5, Pc = 0.4, and Pa — 1. 



At sufficiently high motor susceptibility (5 > for A = 1; Fig. [TOb lower panel), we 
observe a non-monotonic dependence of the tension (— p), i.e. negative pressure, on the 
mean separation i?, clearly indicating the necessity of some kind of phase separation into 
node-rich and node-poor regions. 

For susceptible motors with 5 = 1 (Fig.fTOk). varying the motor activity A affects both the 
short-range (small R) and the long-range (large R) attractions: increasing motor activity (as 
indicated by dashed arrows) leads to stronger localization of individual nodes (upper panel) 
and a larger tension in favor of stronger aggregation (lower panel). In particular at small 
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R where the bonds are buckled, effective attraction arising from motor-driven contractions 
Alogr) dominates, yielding further aggregation of the loosely connected nodes. This 
behavior is most clearly manifested for a sparsely connected network where lack in bond 
constraints allows the formation of isolated aggregates, as shown earlier in Fig. [Tt. 
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FIG. 11. The role of motor kicking rate and effective attraction in aggregation, (a)-(c) Statistical 
measures for the dynamic and structural development at various motor kicking rates. Steady-state 
bond structures (upper) and node configurations (lower) are shown at increasing motor kicking 
rates (d): from left to right n = 0.1,0.2,0.5 and 1 (converted into Tq^/T), and for corresponding 
models without motor-induced short-range attraction at = 1 (e). The remaining simulation 
parameters are Lg — 1.2, /37 = 5, Pc = 0.4, Pa = 1, s — 1, and / — 0.03. 



Using Brownian dynamics simulations, we study the dependence of aggregation strength 
upon motor activity for force-percolating networks {Pc = 0.4). The statistical measures and 
steady-state structures (labeled by Tq^/T) are displayed in Fig.dH Increasing motor kicking 
rate (note A oc /^) apparently enhances the trend toward phase separation (panel d) and 



32 



promotes stronger aggregation (panel c and d), supporting the mean- field prediction. Both 
the fraction of taut bonds (panel a) and the total energy (panel b) decrease with increasing 
kicking rate due to a lower Teff. 

As we pointed out earlier in deriving the effective pair potential, varying motor suscepti- 
bility s affects only the long-range interaction (via Teff), as shown in Fig. [TOb . For a given 
lattice spacing, increasingly susceptible motors (indicated by dashed arrows) drive stronger 
attraction (lower panel) and enhance localization of individual nodes (upper panel), as well as 
yield a stronger tendency for phase separation, as reflected in the increasingly non-monotonic 
dependence of the tension on density change as s increases (lower panel). 




FIG. 12. The stability diagram at various motor susceptibilities. The colored lines represent the 
stability boundaries, solid red for s = 1, dashed grey for 5 = and dotted blue for s = —0.5. Below 
the stability boundaries the pressure exhibits a non-monotonic dependence on particle separation 
indicating the tendency toward phase separation. The instability region (shaded area) extends to 
lower Pc and higher Pa as s increases, suggesting that susceptible motors promote phase separation. 
The remaining simulation parameters are Lg = 1.2, /37 = 5, and A = 1. 



Self-consistent phonon calculations also allow us to determine a quasi- thermodynamic 
stability diagram. In Fig. [12] we present the stability diagram in the Pc-Pa parameter plane 
for susceptible (5 = 1), adamant {s = 0) and load-resisting {s = —0.5) motors. Below the 
stability boundaries (colored lines), the pressure depends non-monotonically on the mean 
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particle separation. This indicates the tendency toward phase separation. Above the bound- 
aries, there are no stable a solutions , or that the tension (or negative pressure) decreases 
monotonically with increasing particle separation. As the motor susceptibility increases, the 
instability region (shaded area) expands toward lower Pc and higher P^, suggesting that 
susceptible motors promote phase separation. Since motor susceptibility affects only the 
long-range interaction via Tgff at high where bond stretching stabilizes finite a solutions, 
small or negative s may lead to an effective repulsion that counteracts the trend of attrac- 
tion thus eliminating the non-monotonicity in pressure, or else destabilizes a solutions by 
offsetting the restoring force. Consequently, the stability boundary at large Pc shifts toward 
lower Pa as s decreases. 

To determine whether the non-monotonicity in pressure indeed corresponds to phase 
separation, we need to examine heterogeneous/site-dependent a solutions; a bimodal dis- 
tribution of stable a values would then indicate that localized dense regions phase separate 
from mobile dilute regions. We hope to investigate this aspect in an upcoming work and thus 
provide more quantitative arguments for the surface tension associated with stretched bonds 
connecting fioppy clumps. Such an analysis should allow us to determine the "droplet" size 
for the condensed phase when bond-constraint-induced arrest occurs. 

To summarize our results in Fig. [13] we outline the diverse patterns formed at various net- 
work connectivity (Pc) for susceptible {s > 0) and load-resisting {s < 0) motors. This figure 
delineates where there are contrasting results from complete MC simulations to those ob- 
tained from BD simulations using the effective pair interaction. As we mentioned before, the 
Monte Carlo schemes explicitly incorporate pairwise anti-correlation by treating move pairs 
as reaction channels, while in the effective Brownian dynamics schemes the total effective 
potential acting on individual nodes from all the neighbors smears out the pair correlation. 
This disparity between the predictions of the two schemes does not affect the active pat- 
terning by susceptible motors (Fig. [T3k). In that case the effective attraction provides the 
dominant mechanism for aggregation and phase separation. MC and BD simulations give 
converging results: Below the percolation threshold (marked by the vertical line) isolated 
fioppy clumps form, whereas at intermediate Pc above the threshold arrested phase separa- 
tion occurs. For load-resisting motors (Fig. [T3b). however, correlation in motion plays a key 
role in structural development, because in this case an effective long-range repulsion dom- 
inates over the short-range attraction and governs the pattern formation. Anti-correlation 
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FIG. 13. Patterns of behavior for susceptible {s > 0) and load-resisting {s < 0) motors at a high 
motor concentration. Typical structures generated by simulations are shown for each situation. 
The horizontal axis indicates increasing network connectivity from left to right. The vertical 
line locates the percolation threshold, (a) For susceptible motors, effective Brownian dynamics 
simulations and Monte Carlo simulations give similar results. At intermediate connectivity above 
the percolation threshold, arrested phase separation occurs, (b) For load-resisting motors, (anti- 
) correlation plays a key role in the active patterning. Macroscopic contraction occurs only in the 
presence of anti-correlation in motion, otherwise connected asters form that cannot collapse. 



in movements gives rise to collapse of the tense intermediates which is essential for active 
contractility /macroscopic contraction.^^ If there is no anti-correlation collapse does not oc- 
cur. As a result, the tense bundles (for non-percolating case) and connected asters (for 
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percolating case) are maintained as long-time steady state structures. At sufficiently high 
connectivity, the bond constraints are too strong to allow significant local force asymme- 
try, thus a balanced tug-of-war between motor-attached filaments prevents the formation of 
heterogeneous cluster structures, and instead leads to a largely homogeneous structure with 
moderate fiuctuations (rightmost regime in both panels a and b), regardless of the motor 
susceptibility or correlation in motion. 

Also we note that the Brownian dynamics simulations for load-resisting motors exhibit 
similar behavior to that coming from uncorrelated isotropic kicks as we studied earlier—: 
Both simulations generate sustained aster patterns that do not collapse, exemplifying a 
negative effective temperature. For susceptible motors, however, the dependence of Teff and 
the effective attraction on the instantaneous local network structure hinders the system 
from achieving global concerted movement. There is thus no spontaneous fiow or oscillating 
behavior that presents for the uncorrelated kick case.— A similar absence of a fiow transition 
has also been found for a contractile nematic model recently studied analytically^^ and 
numerically.— 

IV. CONCLUSION AND DISCUSSION 

We are now in a position to recapitulate how the intricate interplay between local force 
generation, network connectivity and collective action of motors gives rise to regular and 
heterogeneous patterns, arrested coarsening and macroscopic contraction: A sufiicient con- 
nectivity is required for forces to percolate so that local motor-induced stresses and re- 
sulting deformations can propagate significant distances through the network. Given a 
force-percolating structure, downhill-prone motors yield heterogeneous/disordered cluster 
structures, exemplified as an arrested phase separation in the absence of bond or motor rup- 
ture events; whereas load-resisting motors may drive large-scale contraction by surmounting 
a high energy barrier constituted by tense intermediates. Anti-correlation in movements is 
essential for collapse of the intermediate tense states in approach to the eventual large-scale 
contraction; in the absence of correlation in movements, as is the case for Brownian dy- 
namics simulations and for our earlier studies on uncorrelated kicks,— the stretched bundles 
cannot collapse and no contractile structures result, thus the aster pattern is maintained as 
the steady-state structure. 
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The notion of effective interaction provides a natural explanation for the aggregation 
and coalescence of actomyosin condensates: enhanced long-range attraction facilitates ini- 
tial density fluctuations; effective attraction at short distances especially in the buckling 
regime, arising purely from motor-driven contractions, then promotes efficient aggregation 
by drawing nearby nodes or condensates even closer. Whereas an effective repulsion can 
originate from a negative motor susceptibility, and in turn a negative effective temperature 
allows the formation of aster patterns in the absence of correlation in node movements. 

Our finding may suggest a new mechanism for aggregation of active gels: local force 
asymmetry and disorder (structurally inherent or thermally generated) trigger local aggre- 
gation which is further enhanced by an effective attraction due to correlated motor kicks; 
force percolation combined with the tendency to reduce surface tension associated with the 
stretched bonds leads to coarsening of local aggregates; when the balance between local 
bond collapse and neighboring bond stretching is reached, the system forms an arrested 
structure composed of fioppy clumps connected by tense bonds. The pertinent dynamic 
process involves phase separation into node-rich and node-poor regions followed by arrest 
due to bond constraints. 

We also generalize the concept of effective temperature to non-equilibrium many-body 
systems driven by correlated small-step motor kicking events. Explicit simulation tests lend 
support to the validity of picturing such systems as being at an effective equilibrium with 
modified interactions. 
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Science Foundation (Grant PHY-0822283) is gratefully acknowledged. 
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